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We calculate the inclusive heavy quark production cross section for proton-nucleus 
^ | collisions at high energies. We perform calculation in a quasi-classical approximation 

(McLerran-Venugopalan model) neglecting all low-x evolution effects. The derived 
, expression for the differential cross section can be applied for studying the heavy 

> : 

, quark production in the central rapidity region at RHIC. 

(N ■ 

o : 

In this letter we address a problem of heavy quark production in proton-nucleus collisions 
at very high energy. Heavy quarks are a very important tool for studying the properties of the 
strong interactions. At not very high energies the heavy quark mass provides a scale which 
allows one to use the perturbative QCD [1] since the long distance dynamics is effectively 
decoupled Q. However, at high energies which correspond to low values of Bjorken x a 
nucleus becomes a strongly coupled dense partonic system (Color Glass Condensate) with 
large typical transverse momentum Q s determined by density of nuclear partons over the 
nucleus transverse area . Experimental data suggest that Q 2 S ~ 2 GeV 2 for the 

Gold nucleus at x ~ 0.01 7[. It is existence of the strong color field which violates the 
decoupling of the heavv__ quark production subprocess from the dynamics of partons in the 
nucleus wave function 

The Color Glass Condensate starts to play a significant role in scattering processes at 
ow-x since a coherence length of gluons in a nuclear wave function is of the order l/(xM^) 



X 



/y_qu, 







llj which allows them to coherently interact with all nucleons in a nucleus. It was 
argued in refs. ((J 12, ^| that the color field of a nucleus in a low-x regime is given by 



the classical solution to the Yang-Mills equations. It includes all multiple rescatterings of 
a gluons with the color charges of a nucleus 12, l^. However, the quasi-classical approach 



2 



is not sufficient when x < e 1 ' cts . In that case quantum evolution effects become important 
and mUSt be re8UI „ed using the ,ow-, evo.ution e qUati on S SQQQQ. 

In this paper we undertake the first step towards solution of a problem of heavy quark 
production in pA collisions at high energies by deriving the heavy quark production differ- 
ential cross section ()24|) in a quasi-classical approximation which is equivalent to inclusion 
of all multiple rescatterings of a proton with a nucleus. 

The process of heavy quark production at high energies in a quasi-classical approximation 
has three separated in time stages in the nucleus rest frame. Emission of a gluon g by a 
proton's valence quark q v takes much longer time T qv ^ qvg than a subsequent emission of a 
qq pair by a gluon t g ^ q q <C T q v ^q v g- In turn, the time of interaction of a q v gqq system with 
a nucleus is of the order of nuclear length and is negligible as compared to the evolution 
time of the proton wave function T qv -+ qvg ^> t g ^ qq - ^> Ra- Indeed, assume that proton is 
moving in the "+" light cone direction with four momentum p = (j» + ,0,0) and nucleus is 

2 

at rest. Denote the emitted gluon four-momentum by q = {qp + ) ^,p). By uncertainty 
principle emission of a gluon by a valence quark takes time 



T _ = ? = 2,(1 - ,) P+ ^ V (1) 

'' '' q + + q + (jp — q) + + (p — q) — p + — p q 2 q 2 



since the emission of a gluon at high energy is dominated by q <C 1. The Bjorken x is defined 

q 2 

as x = 2q +M N > wnere is a nucleon mass. Therefore 

r„ = ^- » R A ^ r int (2) 

for very low x. The same argument applies to the successive emission of a quark- ant iquark 
pair by a gluon in a proton's wave function: a q v — * q v g fluctuation spans much longer time 
than g — > gg. Therefore, processes in which gluon or heavy quarks are produced in course 
of the rescatterings in a nucleus are suppressed by powers of energy p + . 

Let us choose the A + = light cone gauge. In view of the above argument we can 
separate the process of heavy quark production in nine terms according to the time of gluon 
emission in the amplitude n and in the complex conjugate one T2, time of quark- ant iquark 
emission in the amplitude t\ and in the complex conjugate one £2, and the time of interaction 
which happens at light-cone time T int = x + = 0. In Fig. [T]we show all possible cases. 

To proceed we need to know the light-cone wave functions of a valence quark and of a 
virtual gluon in transverse configuration space. The light-cone wave function of a valence 
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quark in momentum space is given by 

^Ui) = gT a ^ g , (3) 

where q is the gluon's transverse momentum and e A is the gluon's polarization vector. Its 
Fourier image reads 

Averaging square of Eq. (j3J) over quantum numbers of the initial quark and summing over 
quantum numbers of the final quark and gluon we obtain the familiar gluon radiation kernel 
of a dipole model 

where z x and z 2 are the transverse coordinates of the gluon in the amplitude and in the 
complex conjugated amplitude correspondingly. 

Light-cone wave function of a virtual gluon of momentum q reads, see Fig. ^ 

ip g -+ qq -{k,k-q,a) = T, \9~, 9 ($r,r>(k-aq)-e x [r(l-2a) + \]+rS r _ r >m(l+r\)), (6) 

- (k — a q) z + m z - 

where k is the produced quark's transverse momentum, m its mass, a = k + / q + is the fraction 

of the gluon's light-cone momentum it carries, r and r' are the quark and the antiquark 

helicities correspondingly. Eq. © can be written in transverse configuration space using 

modified Bessel functions 

— $((%q ~ Z-i) + a ~~ ^o)) %> — $r,r' — mKi(\x — Xq\ m) [r(l — 2a) + A] 
+-f^o(|^ — £q\ m ) r5 r - r i m (1 + rX) ), (7) 



where x_± and x are the quark's and antiquark's transverse coordinates in the amplitude 
correspondingly, see Fig. and x = Averaging square of Eq. (JZJ) over quantum numbers 
of the initial gluon and summing over quantum numbers of the final quark and antiquark 
[l8 | we find 

fh ( \ ol s 2 ( {x- x^) ■ {y - x^) 

^ g ^qq{z,x,x ,a) = —m -r= ^K l {\x-x \m)K l {\y-XQ\m) 

7T V 37 Xq I | y X_q I 

x[ct 2 + (1 - a) 2 ] + K (\x-x \ m)K (\y- x \ m) V (8) 
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FIG. 1: Diagrams which contribute to the heavy quark production in the light-cone perturbation 
theory. A: r x < 0, t x < 0, r 2 < 0, t 2 < 0, B: t x < 0, r 2 < 0, i x > 0, t 2 > 0, C: r x > 0, i x > 0, 
r 2 > 0, t 2 > 0, D: r x < 0, i x < 0, r 2 < 0, t 2 > 0, E: t\ > 0, t x > 0, r 2 < 0, t 2 > 0, F: r x < 0, 
t x < 0, r 2 > 0, t 2 > 0. Not shown are the complex conjugates D*: r x < 0, t x > 0, r 2 < 0, t 2 < 0, 
E*: t\ < 0, ti > 0, t 2 > 0, t 2 > 0, F*: r x > 0, i x > 0, r 2 < 0, t 2 < 0. Instantaneous interaction of 
a q v gqq system with the nucleus happens at light-cone time Ti nt = 0. The final state is denoted by 
the vertical dashed line at r = oo. 



where y is the quark's transverse coordinate in the complex conjugate amplitude and we do 
not include two delta functions (see (0)) in definition of § g ^ q q- Eq. ® is a special case of 
light-cone wave function of an off-shell gauge boson derived in Refs. 



p is a ; 

3 . y Q. 



Rescatterings of the produced partonic system in a nucleus must be calculated separately 
for each time ordering 20] as shown in Fig. ^ The result can be written in terms of the 
Fourier transformation of the normalized gluon-nucleon cross section jsoL I2I | : 



V(x) 



j2, ^-il-x 1 d(J 



drle 



(9) 



j 



FIG. 2: Diagrams contributing to the time ordering of the diagram F in Fig. ^ 



All diagrams contributing to the time ordering of the diagram F in Fig. are shown in 
Fig. The sum of diagrams a-j in Fig. El is given by 

Cf -v(Q) + t^t(n c - Pi V(x ) 



2N C w 2N r V N, 



= ~ (V(x) - V(0)) + \ (V(x ) - V(0)) - -L (V(x - x ) - V(0)) , (10) 

where x and y are coordinates of the produced quark in the amplitude and in the complex 
conjugate one correspondingly, x is coordinate of antiquark, see Fig.[T] Multiplying expres- 
sion (fTUj) by the nucleus profile function T(6), nuclear density p and the gluon-nucleon cross 
section a |2QJ] we obtain 

- P(x, x ^) = - l -x 2 Ql- l -£ Q 2 s + ^ r2 {x- x,f Ql (11) 



where we follow notations of |22(. The saturation scale Q 2 S in (fTTj) is given by jl 71. \2 

Ql{x) = ^^ P T(b)xG(xA/x 2 ), (12) 
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FIG. 3: Diagrams contributing to the time ordering of the diagram D in Fig. ^ 



where the gluon distribution function in a nucleon reads 

xG{x, 1/x 2 ) 



a s C F 1 
In 



7T 



X 2 /! 2 ' 



(13) 



with /i some infrared cutoff. For spherical nucleus T(b) = 2\J R 2 — b 2 . Assuming that the 
interactions of a proton with individual nucleons are independent we can exponentiate 
the formula (|TT|) to obtain for the diagram F on Fig. 



F = exp {-P(x,x )} = exp < -\xj Q 2 S - \x^Q 2 s + 



{x_ X.Q ) Q 



8N 2 



(14) 



22, 



23 



24] 



This formula coincides with the qqg "propagator" derived in refs. 

Analogously, all diagrams contributing to the time ordering of the diagram D in Fig. ^ 
are shown in Fig. |3] and in Fig. EJ e-j. The sum of diagrams a-i in Fig. |3] and e-j in Fig. |21 
yields 

D = \V(* - *) - ^Vfc) + ivte, - , 2 ) - JL (at. - -1) V(jb ) 



2 K - 2 > 2N r 2N 2 



1 ^-^)-^(o) 



+4^(o) - ^-v(o) - -L (n c - A ) i u n ) 



N 2 



2iV; 



2iV r 



^ (V(x - z 2 ) - 7(0)) + i (Vfe - l2 ) - 1/(0)) - -L 



(V(x-So)-V(O)). (15) 
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Multiplying (JT5Jl by T(6) per and exponentiating we derive 

D = exp{-P(x - z 2 ,x -z 2 )} (16) 



where we used definition ()T4j) . Complex conjugated of F and D can be obtained by replacing 

x^y, and z-y <-> z 2 : 

/J* = exp{-P(y,x )}, (17) 
F* = exp{-P(y-^,x -^)}- (18) 



Diagrams B,C and F, E* have been calculated in |2£j. The only difference is additional 
color factor ~ emerging due to fluctuation of a virtual gluon into quark-anti-quark pair. We 
included this factor in the definition of the wave function (jHJ. We have 

B = e -3fei-&) 2 <^ (19) 
E = e-z£ Q * (20) 
E* = e -*£ Q * (21) 



Finally, it is easy to see that the diagram A in Fig. [T] is equal to 

A = e -\^t3L-y?Ql_ (22) 

Summing up diagrams A - E and their complex conjugates results in the following rescat- 
terings factor 

_l_ e - p fe2o) _|_ e - p (vi2o) — e - p {x-z 2 ,*o- -«a) _ e --P(v-*i>2o-*i) (23) 

Using (jHJ), (jHJ) and (|2HJ) we can write down the inclusive quark production cross section 

da f j2t 12 >2 a s C F z\ ■ z 2 f 2 /■ /• d 2 x d 2 y 
d^ = J dbdZldZ2 — ^[jJ dX0 J da J W~ ^(x-x ,y_-x ,a) 

xe^ ( ^ ) H(x,y,x ,2 1 ,z 2 )(5((x -z-y) + a (x- x )) 5((x - z 2 ) +a(y-x )). (24) 

where h is an impact parameter. This formula is a generalization of result obtained by 
Kopeliovich and Tarasov in Ref. j^jj. 

Formula ()24j) is the main result of our paper. It resums all higher twist effects in the quasi- 
classical approximation which means that we keep all terms proportional to a 2 A 1 / 3 ~ 1 and 



8 



neglect terms suppressed by powers of a s <C 1. We explicitly neglected the low-x quantum 
evolution effects assuming that a s ln(l/a;) <C 1. Therefore formula (}2"4*)) is applicable when 
e -i/a s < X <^1. As resent experimental data on dA collisions at RHIC show, this corresponds 
to the central rapidity region at a/s = 200 GeV j^J. 

Formula (f2"4*|) has been used in Ref. |9J] for numerical calculations of charm production at 
RHIC. It was shown that the charm spectrum obtained according to (}2"4*)) is much harder 
than in naive parton model approach. This is attributed to the presence of a hard 'intrinsic' 
scale Q 2 S . It is clear that the dependence of a heavy quark yield on A is closely related to 
the relation between Q s and m. In the strong color field Q s ^> m the total cross section of 
heavy quark production in pA collisions is proportional to the transverse size of a nucleus 
cr tot ~ A 2 / 3 due to saturation in a nucleus wave function. In the opposite limit Q s <^ m the 
color field of a nucleus is not able to produce heavy quarks from the vacuum in which case 
o"tot ~ A. Therefore, at high energies one expects suppression of the heavy quark yield. In 
the case of charm quark production numerical calculations in [*} show that at y = at RHIC 
the charm quark yield is not suppressed. However at forward rapidities it gets suppressed 
since the nuclear color field strength increases at small x due to quantum evolution. I refer 
the reader interested in phenomenological applications of (|2"4*|) to Ref. for more detailed 
discussion. 

Dynamics of saturated quasi-classical color fields dominates the total multiplicities of 

AA and dA collisions in the central rapidity region at rhic Q,y,y,y,Q- 

The Cronin 

enhancement seen in the data jsjj is produced by multiple rescatterings of a proton [3, [3] 
in a saturated wave function of a nucleus These multiple rescatterings produce 

particle correlations which give a substantial contribution to the elliptic flow phenomenon 
in AA collisions |36[ . Eq. (}2"4*)) can be used for analysis of the heavy quark production in pA 
collisions at the central rapidity region at RHIC. In particular, one can address the question 
of whether formula (J2l| yields the Cronin enhancement of charm production analogously to 
the case of gluon production (sii 

HQ- 
High energy quantum evolution has been neglected throughout this pap er. However, as 
energy /rapidity increases the quantum evolution becomes important [lj, 
gives rise to a number of spectacular effects 



geometric scaling phenomenon 



39 



y 



41 



42 



31 



33l [381 associated with the extended 



Recent results of BRAHMS collaboration 



|43| at RHIC indicate onset of the high energy evolution at rapidities close to the proton 
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fragmentation region in agreement with theoretical predictions. Therefore, generalization of 
(124)1 to include low- a; quantum evolution is an important task which will be pursued in our 
forthcoming publications. 

Finally, we would like to mention an important theoretical question which have not been 
touched in this paper. It is whether eq. (EH) can be reduced to the /cr-factorized form. The 

nnn 

/^-factorization was proved for dilute target regime in 8L 4J, |45J and have been recently 



rederived in a Color Glass Condensate framework in ref 



It was also proved in ref. 



that the inclusive gluon production cross section in pA collisions can be reduced to the 
/cr-factorized form even if the quantum evolution is included. So far all phenomenological 
studies of the heavy quark hadroproduction at high partonic densities q| |3| have been 
based on /^-factorization. Therefore it is important to understand to what extend it can 
be realized at high energies and/or for heavy nuclei. We are going to address this problem 
elsewhere. 
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